Astronomy & Astrophysics manuscript no. HD 144432 


©ESO 2011 


October 18,2011 





HD 144432: A young triple system^ 

A. Mulleiffl, A. CarmonEpE]^ m. E. van den AnckeiSl, R. van Boekeff, Th. Hennin^, and R. Launhardff 

' Max Planck Institute for Astronomy, Konigstuhl 17, D-691 17 Heidelberg, Germany 
e-mail: amueller@mpia.de 

^ ISDC Data Centre for Astrophysics & Geneva Observatory, University of Geneva, chemin d'Ecogia 16, 1290 Versoix, Switzerland 
^ Observatoire de Geneve, University of Geneva, chemin des Maillettes 51, 1290 Sauverny, Switzerland 
e-mail: andres . carmona@unige . ch 

European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching b. Munchen, Germany 
Received date / Accepted date 

ABSTRACT 

We present new imaging and spectroscopic data of the young Herbig star HD 144432 A, which is a well-known binary star with 
a separation of 1.47". High-resolution NIR imaging data obtained with NACO at the VLT reveal that HD 144432 B itself is a 
close binary pair with a separation of 0.1". High-resolution optical spectra, acquired with FEROS at the 2.2m MPG/ESO telescope 
in La Silla, of the primary star and its co-moving companions are used to determine their main stellar parameters, such as effective 
temperature, surface gravity, radial velocity, and projected rotational velocity by fitting synthetic spectra to the observed stellar spectra. 
The two companions, HD 144432 B and HD 144432 C, are identified as low-mass T Tauri stars of spectral type K7V and MIV, 
respectively. From the position in the HRD, the triple system appears to be co-eval with a system age of 6±3 Myr. 

Key words. Stars: pre-main sequence - protoplanetary disks - Stars: atmospheres - Stars: binaries: close - Stars: fundamental param- 
eters - Stars: variables: T Tauri, Herbig Ae/Be 



1. Introduction 



Herbig Ae/Be (HAeBe, iHerbid Il960t iFinkenzeller & Mundj 
11984:) stars are intermediate-mass (2-10 Mq) pre-main sequence 
stars (PMS). They are the higher-mass counterparts of the T 
Tauri stars (TTS), hence fill the parameter space between TTSs 
and high-mass young stars in addressing the question of star for- 
mation as a function of mass. Double or multiple stellar systems 
are of great interest to star formation theories as they provide 
constraints on the components of star formation models, such 
as fragment ation, accretion, N-body dynamics, and orbit mi- 
gration (e.g. lTokovi nin 2008. and references therein). Although 
several imaging survey s have indicated that HAeBe stars have 
close companions (e.g. iLeinert et al ] 1 19971: lPirzkaletal.1119971: 
iKouwenhoven et al] l2005l) . relativelv few studies have con- 
strained the physical properties of these companions. One in- 
teresting question that can be addressed is the age of the HAeBe 
star by constraining the age of its companions. An independent 
and more robust age estimate can be derived from the lower- 
ma ss compan ions . 
' ICorporonI (il99 8) pointed out that the companions of Be stars 
tend to be intermediate-mass stars with spectral types A to F, 
while the companions of Ae stars are low-mass TTSs with spec- 
tral types K to M, which suggests that the mass ratio seems to be 
quant itatively independent of the primary mass (e.g. HogeveenI 
Il992h . 

HD 144432 A is a late A to early F-type star (i The et al.l 

119941) at a distance of 160+^^ pc jvan Leeuw en 2007) that is 
associated with the Sco OB2-2 star forming region. The star is 
surrounded by a protoplanetary disk (e.g. Meeus et al. 2001) and 



* Based on observations made with ESO Telescopes at the La Silla 
Paranal Observatory under programme ID 075.C-0668 and 083.A- 
9013. 



displays signs of active mass accretion. The val ues found in the 
Utera t ure range from 1.8 to 8.5- 1 '^ Mp yr~^ jBlondel & Diid 
120061: iGarcia Lopez et alj 120061: iDonehew & BrittainI 1201 iT 
Frorn__J=3-2_J^CO measurements and line profile fitting. 
Dent et alJ (|2005|) derived a dust mass of 2-10"^ M© (assuming a 
constant gas-dust ratio of 100) and an outer radius of 60+20 AU 
(i - 48° + 10°) for the disk around HD 144432 A where emitting 
gas is present. The system HD 144 432 is known to be a binary 
('Perez et al. 2004; Carmon a et al.l 1200 7). K-hmd observations 
by Perez et al. (2004) revealed a 2.4 m ag fainter stellar com - 
panion at a separatio n of 1.4" (PA ^4°. ICarmona et al]|2007h . 
ICarmona et"an (l2007h reported the detection of Li I (6708 A) 
and Ha emission in the spectrum of HD 144432 B using the 
VLT-FORS2 spec ti-ograph. 

iPerez et alJ (|2004 showed that HD 144432 A and B share the 
same common proper motion and age. In addition, compone nt 
B was ruled out as a background star bv lCarmona et alj (l2007h . 

In this paper, we present NACO observations of HD 144432, 
which resolve HD 144432 B as a close binary star, thereby 
revealing HD 144432 as a hierarchical triple system. Using 
high-resolution optical spectroscopy, we determined the stellar 
parameters of all components of the triple system. 



2. Observations and data reduction 

2.1. NACO 

HD 144 432 was observed i r i the night of July 1 3, 2005, with 
NACO ( Lenzen et all 120031: 'Rousset et all l2003l) . the adaptive 
optics (AO), near-infrared (NIR) camera at ESO's VLT. We ob- 
tained these unpublished observations from the ESO archive. 
The data consist of six imaging observations in three different fil- 
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ters (Ks at /I = 2.18 //m, NB2.12 at ^ = 2.122 yum, and NB2.17 at 
A = 2.166 yum). Each observation consists of eight single expo- 
sures at different sky positions (jitter). Standard data reduction 
(bias subtraction, flat-field correction, bad-pixel interpolation) 
was applied to each frame by employing the available calibration 
frames. The average FWHM of the images is 0.07". A combined 
image was produced by applying the shift-and-add method (a 
two-dimensional cross-correlation routine) to the NB2. 17 obser- 
vation (Fig. [TJ using NACO pipeUne recipes provided by ESCQ. 



2.2. FEROS 

We obtained single-epoch spectra of HD 144432 A (SNR of 350 
at 5500 A) and a combined spectrum of HD 144432 B and C 
(SNR o f 250 at 5500 A) i n the night of June 1, 2009, using 
FEROS ( iKaufer et al.lll999h at the 2.2m MPG/ESO telescope at 
La Silla Observatory in Chile. FEROS covers the optical spec- 
tral range from 3600 A to 9200 A and provides a spectral reso- 
lution of ai48 000. The spectra were obtained using the object- 
calibration mode where one of the two fibers is positioned on the 
target star and the other fiber is fed with the light of a Th Ar+Ne 
calibration lamp. The reduction of the raw data was performed 
using the online data reduction pipeline available at the tele- 
scope. The pipeline performs bias subtraction and flat-fielding, 
traces and extracts the single echelle orders, applies the wave- 
length calibration, and corrects for the barycentric motion. For 
each exposure, it produces 39 individual sub-spectra represent- 
ing the individual echelle orders, as well as one merged spec- 
trum. 

The fiber aperture of FEROS is 2" on the sky. The projected sep- 
aration between the components B and C is 0.1". Therefore, the 
observed spectrum is a superposition of those of the two com- 
panions. In addition, owing to the difference in the brightnesses 
of B and C in the optical with respect to A (5 ~ 5 mag) and a 
seeing at the time of the observations of 0.8", the spectrum of 
the fainter companions is contaminated by the »1.5" apart and 
^20 times brighter primary star. 



3. Data analysis and results 

3.1. Imaging 

To measure the relative position of the companions in the field, 
i.e. their projected distance and position angle with respect to the 
primary star, and its relative fluxes, w e used the IDL-based pro- 
gram Starfinder (iDiolaiti et al.ll2000l) . which is designed to ana- 
lyze AO images using the extracted point spread function from 
the image. The final positional values were derived by averag- 
ing the single values of all individual NACO exposures (48 in 
total). Their errors were derived by computing the standard de- 
viations in the means and taken into account the calibrations of 
IChauvin et aP (l2010h for the plate scale (13.25+0.06 mas) and 
true north orientation (-0.02° + 0. 10°), which were derived from 
a data set taken in August 2005. Table [T]lists the measured dis- 
tances and position angles of the companions. The Kg magni- 
tudes for B and C were obtained from the relative flux measure- 
ments. We computed the zero point of each individual image 
using HD 144432 A as reference (K - 5.888 mag. ICutri et alJ 
120031) and derived 9.09 mag for B and 9.16 mag for C, respec- 
tively. 



Table 1. Relative astrometric measurements of the triple system 
HD 144432, where p is the angular separation, d is the projected 
separation, and PA is the position angle. 



Parameter 


A-B 


A-C 


B-C 


p/["] 


1.465+0.007 


1.468+0.007 


0.102 + 0.001 


d/[A\J] 


234±44 


235+44 


16 + 3 


PA/n 


6.31+0.12 


2.32+ 0.12 


275.90 + 0.28 
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Fig. 1. Near infrared image (logarithmic scaling in brightness) of 
the HD 144432 triple system observed with NACO in July 2005. 
The components are marked by letters in hierarchical order. The 
separation between A and B is 1.47". The components B and C 
are separated by 0.1". 



3.2. Spectral classification of HD 144432 A 

Table|2]lists the derived stellar parameters of HD 144432 A. The 
luminosity, L*, was computed by integrating the total flux under 
the SED. To construct the SED, we used photometric data from 
IdeGeus et al. (1990), Sylvester et al. ( 1996), and Malfait et al| 
(11998 ). and data from the ISO and IRAS point-source catalog. 
We derived an extinction Ay - 0.15 for HD 144432 A by fitting 
a reddened stellar atmosphere model with the T^ff and log g de- 
rived from the spectral fit to the observed UV-optical photomet- 
ric data. The parameters T^ff, log g, and v sin / were computed 
using a self-developed tool for retrieving stellar parameters of 
Herbig stars based on fitting synthetic spectra to the observed 
stellar spectrum. The compu tation of the synthetic spectra was 
carried out using spectrum (Gray & Corbally 1994) with the at- 
LAs9 atmosphere models (Castelli & Kurucz 2004). The grid of 
synthetic spectra has step sizes of 125 K in Teg, 0.15 in logg, 
and 5 kms ' in vsin/. The stellar spectrum was divided into 
four spectral windows (Fig. |2]i. Each window was fitted inde- 
pendently and the individual best fits were averaged to derive 
the final stellar parameters (Table |2]). For a detailed description 
of this procedure, we refer to a forthcoming paper (Miiller et 
al., in preparation). After finding the best-fit value for T^ff and 
log^ at solar metallicity, we tried models with diverse metal- 
licities to find the best fit to the Fel lines. We found that the 
optimal match to the spectrum was provided by solar metallic- 
ity. Figure |2] shows four related windows of the observed stellar 
spectrum of HD 144432 A (black line), which is closely fitted 
by a single star synthetic spectrum (red line). 
We measured the stellar radial velocity, RV, by cross-correlating 
p iljie- Btedi^E s&jiBatt i im with a template spectrum, which was 
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Table 2. Stellar parameters of HD 144432 A, B, and C. 



Component 


A 


B 


C 


Parameter 




Value 




Spectral type 


FOIIIe 


K7V 


MIV 


/ [Lg] 


12.7 


0.50 -0.20 


0.41 -0.16 


Ta/[K] 


7220±115 


4000+250 


3750+250 


logg/ [cms-2] 


3.60±0.15 


4.0±0.5 


4.0+0.5 


/ [Mo] 


1 +0.2 
1.6 -0.1 


0.8+0.2 


0.5+0.2 


R* 1 [Rg] 


2.3±0.5 


1.5+0.7 


1.5+0.7 


V sin i 1 [km s"'] 


75.0+3.5 


55+5 


40+5 


RV llkms.-^] 


-2.8+0.4 


-5+2 


-4+2 


age 1 [Myr] 


9±2 


4 ^^-2 


3^^-i 



1.0 - 
I 0.9 ^\ 



1/ 1rr\ 



4190 4200 . 4210 4220 
Wavelength [A] 



4270 4280 4290 
Wavelength [A] 




4460 4470 4480 
Wavelength [A] 



4560 4570 4580 4590 4600 
Wavelength [Aj 



Fig. 2. The spectral windows used for fitting the synthetic spec- 
trum (red line) to the observed stellar spectrum (black line) of 
HD 144432 A. 



a synthetic spectrum representing the stellar parameters of 
HD 144432 A (Table |2]i. The resulting cross-correlation func- 
tion was fitted by a Gaussian function. The position of the center 
of the Gaussian yields RV. 

From the position of HD 144432 A in the HRD (Fig. [3]l, we 
obtained the stella r mass, rad i us, an d its age using the evo- 
lutionary tracks of ISiess et al.l (2000I). We als o compared with 
the evolutionary tracks of Tognelli et al.l ( 1201 Ih . which yield the 
same results. The derived parameters are in close agreem e nt wit h 
values derived in other publications, e.g.. Dunkin e t al.l (Il997h . 
Ivan Boekel et al.l ( f2005.) . and .Guimaraes et al. (200a). The mea- 
sured effective temperature and surface gravity are comparable 
to those of an FOIIIe main-sequence star 



3.3. Spectral classification ofHD 144432 B and C 

The observed spectrum of HD 144432 B and C is strongly con- 
taminated by HD 144432 A. The contamination varies as a func- 
tion of wavelength, with the blue part of the spectrum being the 
most strongly affected. In addition, the spectra of components 
B and C are superimposed. To derive the spectral type of B and 
C, we compared the continuum-normalized B+C observed spec- 
trum with a continuum-normalized synthetic spectrum made of 
three components: a high-resolution Kurucz theoretical spectra 
of B and C (computed as described in Sec. 13.21 1. and a scaled 
observed FEROS spectrum of HD 144432 A. 

We employed a custom designed interactive IDL softwar^ 
that permitted us to visually compare in real time, after any 
change of parameters, the observed and the synthetic spectra. 




IDL widget software available upon request from A. Carmona 



Fig. 3. The position of the A, B, and C components of 
HD 144432 in an HRD along with the evolutionary tracks of 
ISiess et al ] (l2000h for a metal abundance of Z=0.02. 



The procedure consisted first choosing the spectral types of B 
and C thereby automatically fixing their absolute magnitudes 
(My); the My for the contribution of component A was then var- 
ied to match strength of the features of A; and finally, the v sin / 
for components B and C was set to reproduce the widths of the 
absorption lines. The synthetic comparison spectrum was con- 
structed by scaling the flux of each component by its My, sum- 
ming the scaled fluxes, and normalizing the resultant spectrum 
by dividing it by a second order polynomial fit to the continuum. 
The x^^^ statistic was calculated for each synthetic spectrum. 

We used as a fir st guess a K4V spectrum for B and C 
(ICarmona et al.ll2007h . and searched for a region in the spectrum 
best suited to the spectral classification. We selected the region 
at 7695 - 8100 A because: (/) it has few atmospheric lines; (//) 
it is relatively featureless in the spectrum of A; (///) it has a very 
strong line at 7700 A that allows us to differentiate spectral types 
K and M; (/v) it has several additional weaker absorption lines 
that aid the spectral classification and determination of RV and 
V sin /; and (v) it has a strong absorption line at 7775 A from 
component A that is not blended with absorption lines of K and 
M stars. 

After defining the spectral window for classification, we 
searched for the optimal spectral type pair able to reproduce 
most of the spectral features observed. The best match was given 
by a binary system composed of a K7 V (Tgif =4000 K, log g=4.0) 
and an MIV (reff=3750 K, log^=4.0) star of vsin/=55±5 and 
40+5 kms"', respectively, plus a contribution by A of similar 
magnitude as the K7V component. We checked other regions 
of the spectrum (6020-6420 A, 6330-6740 A, 6655-6760 A, 
7250-7585 A) in addition to the 7695 - 8 100 A region, and con- 
firmed that this binary pair provided a good match (we note that 
the level of A is different for each region). The uncertainty in the 
classification is one spectral sub-class and Alogg=0.5, owing to 
the 250 K and 0.5 grid, used for T^ff and log g. 

Once the spectral types and v sin / of B and C were derived, 
we relaxed the assumption that the RV for both components is 
equal to -3 kms ' (RV of HD 144432 A, Table |2| and shifted 
their spectra until the Gaussian fit to the cross-correlation func- 
tion has its center at km s"' . In this way, we obtained heliocen- 
tric RVs for components B and C of -5 + 2 and -4 + 2 kms"'. 
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Fig. 4. Upper spectra: observed normalized FEROS spectrum at 
the position HD 144432 B and C (in black), and the normalized 
synthetic combined model that most closely describes the ob- 
served spectrum (in red). Lower spectra: Relative contributions 
of HD 144432 B, C, and the contamination from HD 144432 A 
(note the difference in the Y-axis for each spectrum). The verti- 
cal dotted line displays the absorption line of HD 144432 A used 
to scale the flux of HD 144432 A to match the observed spec- 
trum. The strength of the contamination from the A component 
is similar to the flux of HD 144432 B. 

respectively. In Fig. |4] we display the observed FEROS spec- 
trum, the final combined synthetic spectrum of components B, 
C, and A, and the relative flux contribution of each component. 

Finally, we note that to assume that the flux ratio of compo- 
nents B to C is equal to the My ratio, is equivalent to assuming 
that they are at the same distance and have the same extinction. If 
we allow for a different brightness for both components, we find 
that a slightly closer fit to the 7700 A line can be achieved if C 
is slightly brighter (by up to a 0.5 mag). However, it is statisti- 
cally very unlikely that two stars of similar RV, which indicates 
a common membership, and co-eval evolutionary status are so 
close together on the sky, thus it is unlikely that C is brighter 
than B. 

From the measured distances between A, B, and C (Table [TJ 
and their derived masses (Table|2]i, the components B and C orbit 
component A with a period of 2100+600 yr and B and C orbit 
each other with a period of 60+20 yr under the assumption of 
circular orbits and that the triple system is seen pole-on. 

3.4. Signatures of youth in HD 144432 B and C 

The lithium line at 6708 A is visible in the combined spectrum 
of HD 144432 B and C. Because the Hthium line of the B+C ob- 
served spectrum is stronger (EW 0.52+0.02A) than the lithium 
line expected for the combination of a K7V and a Ml V spectra 
(EW 0.42±0.02A), we conclude that HD 144432 B and C are 
likely young stars. 

To constrain the age of the companions from the position 
in the HRD (Fig. O, we estimated their luminosities. We used 
MARCS steflar models (Gustafsson et al. 2008), /Tj-band pho- 
tometry of the NACO images (Sec. 13. II ). and assumed a constant 
extinction of Av=0. 15 (Sec. l3.2l i for all components. We derived 



a luminosity of Lb - 0.50;|;q 2o for the B component and 
Lc - 0.41;^Q L0 for the C component. The derived ages and 
other derived stellar parameters are listed in Table|2] The ages of 
the three stars are comparable taking into account the error bars. 
Therefore, we estimate that the age of the system is 6±3 Myr 
(weighted mean of the individual age values), which is in good 
agreement with the value of 8^j Myr derived bv lCarmona et al.l 
(■2007.) . 



4. Discussion and conclusions 

We have analyzed NIR AO images obtained with NACO of the 
young Herbig star HD 144432 and discovered that it is a triple 
system. The components B and C are located 1 .47" north of the 
primary star, and B and C themselves are separated by 0.1". 
Using high-resolution optical spectra, we have determined the 
stellar parameters Teff , log g, Li,, M*, and v sin / of all three com- 
ponents and derived a spectral type of A9/F0Ve for component A 
and K7V and Ml V for components B and C, r especti vely. This is 
in agreement with the trend suggested by ICorporonl (Il998l) that 
companions of Herbig Ae stars are usually of spectral type K to 
M. 

The derived ages from the position in the HRD of all compo- 
nents are comparable within the error bars and the system age of 
the triple system can be constrained to be 6±3 Myr This implies 
that these stars have a common evolutionary state and therefore 
supports the scenario of formation via cloud fragmentation. 

Acknowledgements. We are very grateful to the referee for comments and sug- 
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